Single-particle and collective excitations in a charged Bose gas at finite temperature 
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The main focus of this work is on the predictions made by the dielectric formalism in regard to the 
relationship between single-particle and collective excitation spectra in a gas of point-like charged 
bosons at finite temperature T below the critical region of Bose-Einstein condensation. Illustrative 
numerical results at weak coupling (r g = 1) are presented within the Random Phase Approximation. 
C , We show that within this approach the single-particle spectrum forms a continuum extending from 

the transverse to the longitudinal plasma mode frequency and leading to a double-peak structure 
as T increases, whereas the density fluctuation spectrum consists of a single broadening peak. We 
, also discuss the momentum distribution and the superfluidity of the gas. 
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The fluid of point-like charged bosons moving in a uniform neutralizing background (CBF or charged boson fluid) is a 
basic model in quantum statistical mechanics 1 , which is thought to be relevant to superconductors in the limit of small 
pairs 2 ^ and to dense plasmas of pressure-ionized helium in astrophysical objects2*£*£. Interest in quantum Bosc fluids 
with long-range interactions has recently extended to Bose-Einstein condensates of atoms with electromagnetically 
^ — 1 induced gravitational-like forces^ and to dipolar Bose-Einstein condensates&2ii£. 

vO ; 

The properties of the CBF at zero temperature have been extensively studied and are quite well known. The 

m : 

strength of the interactions is measured by the dimensionless parameter r s = (47ma 3 3 /3) -1 / 3 , where n is the mean 
density and the Bohr radius, so that the weak-coupling regime corresponds to the limit of high density. This regime 

■ 

has been evaluated in early studies 1 ii2ii2i showing that the ground state is Bose-Einstein-condensed with a quantum 
depletion that increases with r s . At very strong coupling (r s w 160) the CBF undergoes Wigner crystallization 14 , 

o '. 

once the Coulomb repulsions become sufficiently strong to overcome the zero-point kinetic energy associated with 



localization of the bosons on lattice sites. The fluid state at intermediate and strong coupling has been evaluated by 
Quantum Monte Carlo techniques 15 and has been studied by theoretical methods based on the use of Jastrow-Feenberg 
correlated wave functions (see Apaja et a/. 1 - and references therein) or on self-consistent accounts of correlations (see 
Conti et alM> and references given there). 

Relatively less is known on the behavior of the CBF with increasing temperature T towards the critical temperature 
T c for the disappearance of the Bose-Einstein condensate. Hore and Frankel 18 gave analytical results for the longitu- 
dinal dielectric response, and hence for the allowed modes of density oscillation and for the long-range screening of a 
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static point charge, as functions of T within the Random Phase approximation (RPA). This approach is expected to 
be valid at weak coupling and long wavelengths. Hansen et aim have drawn a schematic phase diagram and Bishop 20 
has evaluated the critical behavior near T c within a Hartree-Fock approach using RPA-screened interactions. Davoudi 
et aim have calculated the condensate fraction and the normal and anomalous momentum distributions within a 
Hartree-Fock-Bogoliubov (HFB) approach, adjusted to incorporate the Hugenholtz-Pines relation for the chemical 
potential. 

There are in our view two main conceptual issues that are still worth studying in regard to the physical behavior of 
the CBF at finite T below the critical region near T c . The first concerns the role of the long-wavelength /-sum rules. 
In a normal Coulomb fluid the usual /-sum rule relates the first spectral moment of collective density fluctuations 
(or equivalently the spectral strength of the autocorrelation function of longitudinal current density fluctuations) to 
the total particle number density n and sets the long-wavelength plasma frequency at the value lu p — (47me 2 /to) 1 / 2 , 
independently of temperature 22 . In a superfluid, however, a further sum rule relates the superfluid density n s to 
the spectral strength of the autocorrelation function of the superfluid velocity at long wavelengths 23 . At T = 
one expects that n s — n, but a normal component carrying entropy and viscosity appears and increases in density 
with increasing temperature. In the HFB approach 21 the dispersion relation for quasi-particles was found to tend 
to the value loq = (4nnoe 2 /to) 1//2 at long wavelengths, with no being the condensate density, and this drops with 
increasing temperature up to vanishing at T c . Even though the HFB approach is known to be unreliable in regard to 
conservation principles, it seems proper to further examine the role played at finite temperature by the two distinct 
particle conservation laws in the excitation spectra of a superfluid CBF consisting of a superfluid component and a 
normal component. 

The second issue referred to above concerns the relationship between single-particle and collective density-wave 
excitations in the CBF at finite temperature. A theorem by Gavoret and Nozieres24 states that in the presence of 
the condensate the two types of excitations in a fluid of neutral bosons at T — have the same long-wavelength 
spectrum, apart from a factor. Chiofalo et aim have shown that this theorem also holds for the CBF at T = 0, 
once the Hugenholtz-Pines relation is satisfied. This implies, of course, that the single-particle excitation spectrum at 
long-wavelengths in a Bose-Einstein-condensed CBF at T = has a gap equal to the plasma frequency - in contrast 
to a condensed neutral fluid where the Goldstone theorem requires the presence of a soft acoustic mode related to 
phase fluctuations. The above property of the CBF at T = is again related to sum rules, as discussed by Chiofalo 
et aim*. How this property evolves with increasing temperature still needs clarification. 

In this work we examine the answers given to these questions by the so-called dielectric formalisni2Si2L2&, which 
was earlier developed for a neutral Bose gas and should be eminently suitable for a fluid with long-range interactions. 
This microscopic approach is based on a perturbative expansion and treats on the same footing the single-particle 
propagators and the density fluctuation propagator. The plan of the paper is briefly as follows. In Section 2 we 
recall the essential results of the perturbative expansion for the structure of the propagators, outline their relation to 
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the linear response functions and to the momentum distribution of the fluid, and discuss poles and zeroes of these 
functions in the complex frequency plane. Section 3 reports illustrative numerical results that we have obtained within 
the RPA for the CBF at weak coupling (r s = 1) and discusses the RPA predictions on its superfluidity and momentum 
distribution. Section 4 concludes the paper with a summary of its main results and a critique of the RPA. 

2. DIELECTRIC FORMALISM FOR A CHARGED BOSE GAS 

We consider a three-dimensional homogeneous fluid of spinless Bose particles interacting through Coulomb forces 
described in Fourier transform by the potential Vk = 4-7re 2 /fc 2 for k ^ and Vk — for k = 0. The latter condition 
takes care of charge neutrality as ensured by the presence of a uniform background at the same density n. 

The fluid is treated at finite temperature T by the imaginary-time technique. The collective density-fluctuation 
correlations are described by the two-particle Greens function 

X (fc,T) = -(T T p k (T)p_ k (0)). (1) 

The density fluctuation operator ,o k is given as p k = X) q a q- k /2 a q+ k / 2 m t erms °f t ne creation and destruction 
operators a k and a k for particles with momentum k, which satisfy the usual bosonic commutation relations. The 
presence of a Bose-Einstein condensate made of Nq particles induces hybridization between x(fc,r) and the single- 
particles Green's functions G^(fc,r), which are defined as 

Gi(k,r) = -(T T bi(rKH0)). (2) 

Here = & k if a = +1 and = &^L k if a = — 1, with 6 k = a k — v^o^k- 

We report below the main results of the diagrammatic analysis for the structure of the Greens functions in Fourier 
transform with respect to the imaginary-time variable, following the work of Szepfalusy and Kondor— . The Greens 
function G^(fc,T) and its Fourier transform G^,(k,iuj n ) are related by 

Gi(k,r) = kBTj2<*P(iUnT)G a (k,iu; n ) (3) 

(and similarly for x), where uj n — 2nnkBT are the Matsubara frequencies (h = 1). 

2.1. Results of the perturbation expansion 

In the presence of the condensate every dynamical quantity can be separated into a "singular" reducible part 
(denoted by an index s) and a "regular" irreducible part (denoted by an index r). One can also extract the proper 
part of each quantity (indicated by a superscript ~) with respect to interaction lines. In particular, the density-density 
response function takes the form \ = ~ v kx) where \ = + X • Thus the dielectric function e, defined as 
usual through the relation 1/e = 1 + VkX [22], is determined by the sum of two polarizability functions, 

e = I - v k x {s) - v k x {r) . (4) 



In this formalism the CBF is a two-component fluid, in which the singular component could be viewed as the condensate 
and the regular one as the noncondensate2& (see the RPA results in Section 3 below for an explicit approximate 
realization of this two-fluid model, in which the condensate density will turn out to be the same as the superfluid 
density). One can also introduce at this point a dielectric function for the regular component, = 1 — VkX^ ■ 
The structure of the single-particle Greens functions is similarly obtained with the help of the Dyson equations 
introducing the self-energies, Eg. One finds 



G* 3 = - 
D 



(•5) 

where A Q are the vertex functions, N@ = S^(aiuj n + e k ) with e k = k 2 /(2m) referred to the chemical potential ^t, and 

D = (iw n - e k - E+)(iw„ + e k + EI) + (6) 

The proper self-energies E^ are introduced by setting E^ = Eg + A a v k A^ / e^ r \ and the proper single-particle Green's 
functions are given by & a = N^/D where D is given by Eq. (6) with E^ replaced by Eg\ Finally, the vertex functions 
relate the proper singular susceptibility to the proper self-energies through x^ = A a G^A^ (here and in the 
following summation over repeated indices is implicit). 

The presence of the condensate thus induces a close relationship between charge density fluctuations and single- 
particle excitations. Using the equality 

De = De^ (7) 
the expressions for \ an d can be written in the form 



*M A a NfA? 



and 



G p a = — + af3 ^ , (9) 

where 

e = e W + v k A a N^/D (10) 

and 

D = D + v k A a N^/e ( - r l (11) 

These expressions will be used later below to discuss the poles and zeroes of the various response functions. It may 
also be explicitly remarked that the vertex functions vanish at T c , so that the single-particle and density- fluctuation 
spectra become uncoupled for T >T C . 
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2.2. Superfluidity and momentum distribution 

We pause at this point to set out the expressions for some further physical properties of the Bose fluid that can 
be obtained from knowledge of its Green's functions. The autocorrelation function Xv,v s of the superfluid velocity v s 
is obtained from the component of the matrix of single-particle Green's functions on the Pauli matrix (the 
2x2 identity matrix), 

k 2 



Xv s V s o 

Ul Tig 



G (0) . (12) 



This is related to the autocorrelation function x<j>4> of phase fluctuations as Xv s v s = k X^/vn ■ The superfluid mass 
density p s — TJiTig follows from \v s v s by means of the sum rule-* 



— = lim Xv„v s (k, 0) = lim / — 



50 da; ImXv s v s (k, uj) 



(13) 



The last identity in Eq. (13) emphasizes the meaning of this relation as an /-sum rule on the superfluid component at 
each temperature T, on account of the expression J s = n s v s for the superfluid current density J s and by comparison 
with the /-sum rule on the total density n, 

dijj lmxjj{k,uj) 

Here xjj is the longitudinal current-current response function. Of course, the continuity equation for density fluctu- 
ations relates xjj to the density-density response function \- 

Equation (13) leads to the Josephson relation for the superfluid density, 

"/c7(2m) -/! + £+ -£+ 



— = lim 

no k^o 



(15) 



k 2 /(2m) 

An alternative definition of the superfluid mass density is obtained 29 by writing it as p s = ran — p n , where the mass 
density p n of the normal component is determined by the autocorrelation function of the transverse current density 
Xt as 

p n = lim xr{k,0). (16) 

k— >0 

As discussed in the book of Forster— , a proof of the equivalence between the definition of superfluid density that 
is given in Eq. (15) from properties of the single-particle excitation spectrum and the definition that follows from 
Eq. (16) in terms of the transverse current correlations is still a question of continued interest. 

Turning to the one-body momentum distribution function n(k) of the Bose fluid, this is related to the single-particle 
Green's function by 

n(k) = -fcflTy^ exp(iw n T))GX(k, iui n ) (17) 



where 7/ is a positive infinitesimal. The summation can be written in the form of a contour integral by means of the 
relation 

E<«p<'<«^:<*.-") = sib jf ^Xl-i cilt '"' (18) 

as is depicted in Fig. 1. Equation (18) follows from the fact that the poles of [exp(z/fcsT) — 1] _1 lie at the points 
iu) n on the imaginary axis in Fig. 1. The contour of integration can be conveniently deformed from c\ to C2 under 
conditions such that G+(fc, z) only has single poles with a small imaginary part (see Fig. 1 again). In particular this 
is true at T = for the CBF at long wavelengths, where the poles lie at the plasma frequency uj p — (Anne 2 /ro) 1 / 2 
and infinitesimally away from the real axis 2 ^. In this case we find 

v k A 2 



lim n(k) — 



(19) 

fe^O 



and comparison with the exact result 



obtained by Chiofalo et alJ* yields 



mn ujp 

te n( * ) = -2*5- (20) 



lim A 2 (fc,Wp)/eM(fc,Wp) = n (21) 

at T = 0. A similar approach can be used to infer that in a neutral Bose gas the momentum distribution scales as 
k~ x at low momenta. 

2.3. Spectral structures 

We return to Eq. (7) and notice that, if the dielectric function e of the fluid vanishes at some frequency (z±, say) 
in the complex frequency plane, the function D vanishes at the same frequency. Indeed, from Eq. (10) e^ r ' will not 
vanish there as long as the vertex functions are different from zero. A "ghost" of the collective oscillations of the 
particle density, which are determined by the zeroes of a 2 ^, will therefore be present in the single-particle excitation 
spectrum below T c . In particular, in the long wavelength limit at T = both spectra are delta-function peaks lying 
at the plasma frequency uj v and the Gavoret-Nozieres theorem holds, as demonstrated by Chiofalo et a/. 25 . The 
hybridization of the poles in \ an d m the G's can be understood by viewing the single-particle spectrum as associated 
to excitation of a particle out of the condensate, which should be driven to oscillate by the Coulomb forces and in 
turn drive an oscillation of the whole fluid at the plasma frequency. 

However, the two types of spectra should be expected to become substantially different at T > 0, when a normal 
(non-superfluid) component of the fluid emerges and permits different damping mechanisms for density fluctuations 
and for single-particle excitations through the dissipation of momentum and energy. Damping by excitation of 
longitudinal currents is operative for the former, and by excitation of both longitudinal and transverse currents for 



the latter. Transverse current fluctuations correspond to a pole in the dielectric function^ at some frequency (22, say) 
in the complex frequency plane. From Eqs. (7) and (10) this can correspond to a zero in D and/or to a pole in at 
the same frequency. One may thus envisage the spectrum of single-particle excitations at finite temperature as forming 
a continuum that extends from a lower edge determined by the onset of transverse currents and is surmounted in its 
upper part by a ghost due to the coupling to collective density fluctuations. The latter, on the other hand, correspond 
to longitudinal current fluctuations and are not coupled to transverse current fluctuations. We shall explicitly see in 
Section 3 below how these expectations are confirmed by the RPA calculations. 

One last point to be made concerns the dielectric function e^ r \ which gives the dielectric response of the "regular" 
component of the fluid. From Eq. (11) we see that a zero of this function at some frequency (23, say) in the complex 
frequency plane corresponds to a pole in D at the same frequency. It is tempting to view this fact as expressing the 
effect of transverse single-particle motions in driving the regular component. It should also be noted that the form 
of e^(k,z) may depend on the collisionality of the gas, which may induce a dependence on kl where I is the mean 
free path between collisions (see the discussion given by Pippard^I for an analog in the gas of conduction electrons 
in metals). This may even affect the ability of the regular component to screen the density fluctuations and the 
single-particle motions. 

We proceed from this heuristic overview to examine how the expected spectral features are brought out in the 
approximate theoretical model afforded by the RPA. 

3. THE RANDOM PHASE APPROXIMATION 

As developed by Hore and Frankeli& and by Szepfalusy and Kondor— , the RPA for a Bose fluid at temperature 
T < T c adopts the ideal Bose gas as reference and accounts for the interactions by a mean-field treatment neglecting 
all short-range correlations. In a fluid of charged particles these approximations are expected to be valid at long 
wavelengths (k — > 0) and for low values of the Coulomb coupling strength r&&. 

As a first step, therefore, the chemical potential is set to zero and the proper susceptibilities are taken to be those 
of the ideal Bose gas, that is 

X« M = ^% (22) 

for the condensate and 

r (r)/> \ /(q) - /(q + k ) 

X () (k,u)= ^ , -T (23) 

for the thermal cloud, with /(q) = [exp(e g /fcsT) — 1] _1 being the statistical distribution function of the ideal Bose 
gas below T c (we are for simplicity omitting to write an infinitesimal imaginary part in the denominators). The 
function e^ r \k,uj) accordingly describes the dielectric response of the non-interacting thermal cloud and keeps no 
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account of interparticle collisions: in the language of quantum hydrodynamics, the thermal cloud is taken to be in 
the collisionless regime. 

At the next step one sets A Q = A Q = ^/uq and = 0. One accordingly finds 

D(k,u) = lo 2 - e\ (24) 

and 



The approximations listed above finally yield 



e(k, U )=eM(k, U )-?P%, (26) 

LO A - El 



D(k, lo) = lo 



2 _ 2n VkSk 



and 



G+(fc,w) = 



lo + e k + nQVk/e^ik,^) 



(27) 



(28) 



D{k,w) 

In the RPA the self-energies are determined by Coulomb interactions with the condensate screened by the ideal 
thermal cloud and the transverse current modes lie at frequency Ek- 

The Bogoliubov approximation is recovered by setting e^ r \k,Lo) — 1. As discussed by Szepfalusy and Kondor— , 
such a Bogoliubov regime can be met in a neutral Bose gas under somewhat special conditions at intermediate 
temperatures. 



3.1. Excitation spectra 



It is immediately seen from Eqs. (26) and (27) that in the RPA the single-particle Green's functions and the 
density-density response function have a pole at the same frequency, which is determined by the zero of e[k, z) in the 
complex frequency plane. In addition the single-particle Greens functions may have a second pole determined by a 
pole in e^ r \k, z). 

The calculation of the susceptibility of the ideal Bose gas^S yields 



and 



Rex (r) (k,z) 



4ir 2 k 



dqqn(q) In 



m 2 z 2 - (qk - ifc 2 ) 2 



2 z 2 -(gfc + ifc 2 ) 2 



Imx^ r '(k, z) — — r In 
y ' 47rfc 



l-exp[-(z-£ fc ) 2 /(4fc B r £fe )] 



l-exp[-(z + £fe ) 2 /(4fc B T£ fe )] 



(29) 



(30) 



In the limit fcA^s ^S> 1, where XdB = (l/2mfcsT) 1 / 2 is the thermal de Broglie wavelength, £( r ) acquires the single-pole 
structure 

* w (M)^ 2( " a -"° 2 )efc . (3D 

Under conditions such that the imaginary part of x can t> e neglected, the zero of e(k,uj) yields the dispersion 
relation 

E k = y/el + 2nov k e k /e^(k,E k ). (32) 

and in particular, in the extreme low-temperature limit where Eq. (31) holds at all wave numbers, one finds E k — > 
\j2nv k e k + - ► w p at long wavelengths. 

The dispersion relation obtained from Eq. (32) in the RPA for r s = 1 is shown in Fig. 2 at increasing values of the 
temperature (in units of K 2 / {kBma 2 B ), such that T c = 1.27). A fc 2 term is present at finite T, where one findsi^ 

El =lu 2 p + A{n-n )k B Ts k + el (33) 

with a numerical coefficient A = 3£(f )/C(l) given in terms of Riemann zeta functions. The effect of short-range 
correlations beyond the RPA is to introduce a term of order fc 2 even at T = with a negative coefficient, which is 
responsible for the presence of a roton minimum in the dispersion relation. 

Figure 3 shows the spectrum of collective density oscillations at a relatively small value of the wave number 
(k = I/as), as calculated in the RPA for the weakly coupled CBF (r s = 1) at increasing values of temperature. The 
spectrum of single-particle excitations is shown in Fig. 4 at the same wave number and temperatures. At such wave 
number the plasmon excitation is still seen in Fig. 3 as a sharp low-temperature peak, which with increasing T moves 
upwards as predicted by Eq. (33) and at the same time broadens. Figure 4 shows that the spectrum of single-particle 
excitations is very substantially different: its main peak is definitely broader at the lowest temperature, and as T 
rises this peak broadens much more rapidly and at the same time acquires a low-frequency plateau which at high 
temperature is developing into a second peak. The sharp feature at the bottom edge of the spectrum lies for all 
values of T exactly at the frequency e k of transverse current excitations, which in the units used in Fig. 4 is e k — 0.5. 
The shift of oscillator strength from the higher to the lower peak in the single-particle spectrum with increasing T 
signals the progressive reduction of the condensate and hence of the hybridization with collective density fluctuations 
hybridization vanishes with the condensate at T c . 

3.2. Superfluidity 



From Eq. (12) we find that the autocorrelation function of the superfluid velocity within the RPA takes the form 
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The sum rule in Eq. (13) (or equivalently the Josephson relation in Eq. (15)) then yields 

p s = mno, (35) 

that is, within the RPA the superfluid density and the condensate density coincide. 

The same result is found from the second definition of the superfluid density based on Eq. (16) for the mass density 
of the normal component. From the dielectric formalism we have \t = XT and within the RPA the latter function is 
approximated by the transverse current response of the ideal Bose gas: 

The imaginary part of the expression in Eq. (36) can be evaluated analytically by the method used for the ideal Fermi 
gas by Nifosi et alM. We set lm.XT(k, u>) = A^{uj) — Ak(— u>), where 

\2 



q 

By direct integration we obtain 



2 



AkH = ^-^^PolyLog 

Z7T 



2, exp 



(w - Efc) 



2 



(38) 



where PolyLog[n, z] denotes the n th -order polylogarithm function of the variable z. Numerical integration yields 
Xt(^) 0) and hence with Eq. (16) we find p n = m(n — no) m the RPA. These results extend to finite temperature the 
findings of Alexandrov and Beere^, who have calculated the transverse response using the Bogoliubov theory. 

3.3. Momentum distribution 

The momentum distribution in the RPA is obtained from Eqs. (17) and (18) using the result in Eq. (28) for the 
Green's function G\. It is useful in performing the integration in the complex frequency plane to rewrite Eq. (28) in 
the form 



G +(k z) - 9+{k,z) _ gJJ^z) 



where f(k, z) — \J z 2 + 2nov k e k [l — l/e( r )(fc, z)], E^ = ^/2noe k v(k) + e\ is the Bogoliubov dispersion relation, and 

, e k +n v k /e^(k,z) z 
5±(M) = 2Ef ± ^M)- (40) 

The integration in Eq. (18) can be carried out by deforming the contour from c\ to ci as depicted in Fig. 1, provided 
that the spectral broadening due to thermal excitations can be neglected. The result is 



n(k) 1 



\df/dz\ z=l 



g+(k,E k ) g-{k,E k ) 



exp{E k /k B T) - 1 1 - cxp{-E k /k B T) 



(41) 



where E k is the dispersion relation given in Eq. (33). A plot of Eq. (41) is shown in Fig. 5. 
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It is seen from Fig. 5 that Eq. (41) misses the subleading term in the momentum distribution at low momenta, 
which is known to be proportional to 1/fc with a positive coefficient^ and is present in the HFB approach—. It may 
also be recalled that the HFB calculations yield a sizable quantum depletion of the condensate in the CBF already at 
r s = 1, as well as a sizable upward shift of the critical temperature over the ideal-gas value 21 . In principle these effects 
could be examined in an extended RPA approach, where a measure of self-consistency may be included through a 
redetermination of the condensate fraction from the one-body momentum distribution function. Such an effort seems, 
however, hardly worthwhile. 

4. CONCLUDING REMARKS 

In conclusion, we have studied the dynamical properties of a weakly coupled charged Bose gas at finite temperature 
by means of the dielectric formalism. The formalism allows one to treat at the same time and on the same footing 
both the collective density fluctuations and the single-particle motions, keeping account of the hybridization between 
one-body and two-body Green's functions typical of a Bose-Einstein-condensed fluid. The formalism could also be 
extended to treat spinorial Bose gases. 

We have illustrated the dielectric formalism by numerical calculations carried out within the random phase approxi- 
mation, in which the dynamical effects of the noncondensate enter through the screening of the Coulomb interactions. 
We have seen that this approximation yields a qualitatively correct view of the spectra of collective and single particle 
motions, in which the latter at finite temperature are directly coupled to both the longitudinal and the transverse 
current density fluctuations. The spectral strength of the single-particle motions gradually shifts from the former to 
the latter as T is raised towards T c . It would indeed be absurd that, as T is near to entering the critical region, the 
excitation of a single particle in the presence of a disappearing condensate should still be expected to require that the 
whole gaseous cloud be set into oscillation at the plasma frequency. 

The main limitation of the dielectric approach is that the quantum fluid is treated as composed of a condensate 
and a noncondensate, rather than as consisting of a superfluid and a normal fluid, so that the superfluid fraction is 
an outcome of the calculation rather than its basic input. This viewpoint is, of course, intrinsic to the Bogoliubov- 
Hugenholtz-Pines theoretical framework. In the RPA the superfluid fraction turns out to be the same as the condensate 
fraction, and the shortcomings of the one-body momentum distribution that it predicts have been pointed out in the 
body of the paper. An account of this and other quantitative defects of the RPA predictions, such as the lack of a roton 
minimum in the dispersion curve of collective excitations, would require a refined treatment of the proper self-energies 
and the inclusion of short-range correlations between density fluctuations. The so-called one-loop approximation 2 ^ 
could be useful in these respects. Further spectral structures may also be expected to arise— 
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FIG. 1: Integration contours in the complex frequency plane. 
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FIG. 2: RPA dispersion relation Et (in units of h /(ma B )) as a function of wave number k (in units of 1/as) for the CBF 
at r s — 1, at the values of the temperature T indicated in the figure (in units of h 2 /{kBma 2 B ) - in these units the critical 
temperature for Bose- Einstein condensation of the ideal Bose gas is T c = 1.27). 
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FIG. 3: Spectrum of collective density excitations at kas = 1 for the CBF at r s = 1 in the RPA, at the values of the temperature 
T indicated in the figure (in units of h 2 /(kBma 2 B )). 
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FIG. 4: Spectrum of single particle excitations at fcas = 1 for the CBF at r s = 1 in the RPA, at the values of the temperature 
T indicated in the figure (in units of h 2 / '(fcsmag)). 
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FIG. 5: Momentum distribution k 2 n(k) as a function of wave number k (in units of 1/as) for the CBF at r s = 1 in the 
RPA, at the values of the temperature T indicated in the figure (in units of h 2 / {kBma 2 B )). The classical Maxwell-Boltzmann 
distribution is also shown. 



